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Abstract. In this paper we study heat kernels associated to a Carnot group G, endowed with a family 
of collapsing left- invariant Riemannian metrics cr^ which converge in the Gromov-Hausdorff sense to a sub- 
Riemannian structure on G as e — > 0. The main new contribution are Gaussian-type bounds on the heat 
kernel for the metrics which are stable as e — > and extend the previous time-independent estimates in 
1161 . As an application we study well posedness of the total variation flow of graph surfaces over a bounded 
domain in (G, Cj). We establish interior and boundary gradient estimates, and develop a Schauder theory 
which are stable as e — >■ 0. As a consequence we obtain long time existence of smooth solutions of the 
sub-Ricmannian flow (e = 0), which in turn yield sub-Ricmannian minimal surfaces as t — ^ oo. 



1. Introduction 
Models of image processing based on total variation flow 

dtu = (iw(Vu/|Vu|) 

have been first introduced in by Rudin, Osher, and Fatemi in |43| . in order to perform edge preserving denois- 
ing. We refer the reader to the review paper |19| where many recent applications of total variation equation 
in the Euclidean setting are presented. The flow t — >■ u{-,t) represents the gradient descent associated to 
the total variation energy / |Vm| and as such has the property that both the total variation of the solution 
t — > u{-,t) and the perimeter measure of fixed level sets = const} are non-increasing in time. Aside 

from its usefulness in image processing, the flow also arises in connection with the limit of solutions of the 
parabolic Laplacian dfUp = div{\S/up\P^'^'S/up) as the parameter p — >■ 1+. In the case where the evolution 
of graphs St — {{x,u{x,t))} is considered, i.e. dtu — div{S/u/ ^1 + |Vup), then in both the total variation 
flow and the closely related mean curvature flow dtu = \/l + \Vu\'^div(Vu/ y^l + |Vup), given appropriate 
boundary /initial conditions, global in time solutions asymptotically converge to minimal graphs. For further 
results, references and applications of the Euclidean total variation flow we refer the reader to [T], [2], [5] 
and for an anisotropic version of the flow to |36| . 

In this paper we study long time existence of graph solutions of the total variation flow in a special class 
of degenerate Riemannian ambient spaces: The so-called (sub- Riemannian) Carnot groups [26], |46| . Such 
spaces are nilpotent Lie groups endowed with a metric structure (G, ctq) that arises as limit of collapsing 
left-invariant tame Riemannian structures (G, Ce). Our motivations are twofold: 

(a) On the one hand, minimal surfaces and mean curvature flow in the Carnot group setting (and in 
particular in the special case of the Heisenberg group) have been the subject of numerous papers in recent 
years, leading to partial solutions of long-standing problems such as the Pansu conjecturqj. Despite such 
advances little is known about explicit construction of minimal surfaces. Since the asymptotic limit of 
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graphical mean curvature flow provide minimal surfaces in the Riemannian setting it is only natural to 
follow the same approach in the sub-Riemannian setting. However this avenue runs into considerable (so 
far unsolved) technical difficulties due to the combined effect of the degeneracy of the metric and the non- 
divcrgcncc form aspect of the relevant PDE. In the present work wc provide the basis for the construction 
of sub-Riemannian minimal graphs through the asymptotic behavior as t — > oo of solutions of the sub- 
Riemannian total variation flow. For previous work on a similar theme see the work of Pauls |40j and Cheng, 
Hwang and Yang [14]. 

(b) A new class of image processing models based on the functionality of the visual cortex have been 
recently introduced in Lie groups. With a generalization of the classical Bargmann transform studied by 
Folland (see [27]) which lifts a function to a new one, deflned on the phase space, an image can be lifted to 
a Lie group with a sub-Riemannian metric. The choice of the Lie group depends on the geometric property of 
the image to be studied: in [18] , and [23] , (JJ 2D images are lifted to surfaces in the Lie group of Euclidean 
motions of the plane to study geometric properties of their level lines, in |25| cardiac images are lifted in 
the Hciscnbcrg group to study their deformations, in [3] moving images are lifted in the Galilei group. The 
image processing is performed in these groups with algorithms expressed in terms of second order subclliptic 
differential equations. In particular sub-Ricmannian mean curvature flows or total variation flows can be 
applied to perform image completion or impainting in this setting. 

Our approach to the existence of global (in time) smooth solutions is based on a Riemannian approximation 
scheme: We study graph solutions of the total variation flow in the Riemannian spaces (G,CTe) where G is a 
Carnot group and is a family of tame Riemannian metrics that 'collapse' as e — >■ to a sub-Riemannian 
metric ctq in G. The main technical novelties in the paper are a series of different a-priori estimates, which 
are stable as e — J- 0: 

(1) Heat kernel estimates for sub-Laplacians and their elliptic regularizations (see Proposition 12. 2[) . 
These estimates are established in the setting of any Carnot group and provide a parabolic counterpart 
to the time- independent estimates proved by two of the authors in [16]. 

(2) Uniform Schauder estimates for second order, non-divergence form subelliptic PDE and their elliptic 
regularizations (see Proposition 14. 4[) . These estimates are established in the setting of any Carnot 
group. 

(3) Interior gradient estimates for solutions of the total variation flow (sec Proposition 13. 2[) . These 
estimates are established in the setting of any Carnot group. 

(4) Boundary gradient estimates for solutions of the total variation flow (sec Proposition 13. 4p . These 
estimates are established only in the setting of Carnot group of step two. 

We remark that (3) and (4) can also be proved, with similar arguments, for solutions of the sub-elliptic mean 
curvature flo-vj^. The limitation to Carnot groups of step two in (4) is due to the fact that for step three 
and higher there no known suitable barrier functions, or in other words, the class of known explicit minimal 
graphs is not sufficiently large. 

In order to state our results we need to introduce some notation. 

1.1. Carnot group structure. Let G be an analytic and simply connected Lie group with topological 
dimension n and such that its Lie algebra Q admits a stratification Q = © ... © y, where [y^, V^] = 
yj+i, if j = 1, . . . , r- 1, and [V'', y] = 0, A: = 1, . . . , r. Such groups are called in [26], [28], and [46] stratified 
nilpotent Lie groups. Set H{Q) = V^, and for any x G G we let H{x) = xH(Q) = span(Xi, . . . , Xm){x). The 
distribution x — > H{x) is called the horizontal sub-bundle H . 

The metric structure is given by assuming that one has a left invariant positive definite form ctq defined 
in the sub-bundle H. We fix a orthonormal horizontal frame Xi, . . . ,Xm which we complete it to a basis 



•^For more references on the sub-Riemannian mean curvature flow see [8], |22| , |f 2| and |35 | 
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(Xi, . . . , Xn) of Q by choosing for every fc = 2, . . . , r a basis of 14 . If Xi belongs to Vk, then we will define 

its homogeneous degree as 

(1.1) d{i) = k. 

Wc will denote by xX = ^iXi a generic element of Q. Since the exponential map exp : — > G is a 

global diffeomorphism we use exponential coordinates in G, and denote x — {xi , . . . , 2;„) the point exp {xX^ . 
We also set xh = (xi, . . . ,Xm)- Define non-isotropic dilations as Ss{x) = (s'^^^^a;^), for s > 0. 

We let Vo = {Xi, . . . , Xm) denote the horizontal gradient operator. If ^ e C°°{G) we set Vo(/> = 
E"iX,</>X, and |Vo</>P = J:T=liX^<P?■ 

We denote by {Xi, . . . ,Xn) (resp. (X[, . . . ,Xj-J) the left invariant (resp. right invariant) translation of 
the frames (Xi, . . . , X„). 

The vectors Xi,...,Xm and their commutators span all the Lie algebra G, and consequently verify 
Hormander's finite rank condition ([35]). This allows to use the results from and define a control 
distance d^^x^ y) associated to the distribution Xi , . . . , Xm, which is called the Carnot-Carathcodory metric 
(denote by d^.o the corresponding right invariant distance). We call the couple (G, do) a Carnot Group. 

We define a family of left invariant Riemannian metrics , e > in C/ by requesting that 

{XI, . . . , X^} := {Xi, . . . , Xm, eX,n+i, . . . , eXn} 

is an orthonormal frame. We will denote by d^ the corresponding distance functions. Correspondingly we 
use Vc, (resp. V^) to denote the left (resp. right) invariant gradients. In particular, if g G°°(G) we set 
= E" 1 and |V,</>|2 = S^^^(X|<^)2. 

We conclude by recalling the expression of the left invariant vector fields in exponential coordinates (see 

m) 

r 

(1-2) X, = d,+ J2 E Pik(^)d,, 

k=d{i) + l d{j)=k 

where pli^{x) is an homogeneous polynomial of degree k — d{i) and depends only on Xh, with d{l) < d{h) < 
k - d{i). 

1.2. The total variation. The total variation flow is characterized by the fact that each point of the 
evolving surface graph moves in the direction of the upward unit normal with speed equal to the mean 
curvature times the volume element. In the setting of the approximating Riemannian metrics (G,CTe) and in 
terms of the functions t — u{-,t) : £7 C G — M describing the evolving graphs, the relevant equation reads: 

(1-3) ^ = - E^Kitr) ^ ^ 4(v.-e)x.^x>, 

for X G O and t > 0, with Ue(x,0) ~ is the mean curvature of the graph of Ue{-^t) and 



W^^l + \y,u,\^ = 1 + ^(X|7.,)2 and al^iO ^ W^^ % - - 



i = l 

for all ^ G M". In the sub- Riemannian limit e = the equation reads 

i—l ^' 

for X G and t > 0, with 0) — ^{x)- 



v/l + |Vou| 
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We will be concerned with uniform (in the parameter e as e — > 0) estimates and with the asymptotic 
behavior of solutions to the initial value problem for the mean curvature motion of graphs over bounded 
domains of a Carnot group G, 



(1.5) 




in Q = X (0, T) 
on dpQ. 



Here dpQ = (Q, x {t = 0}) U [dO. x (0,T)) denotes the parabolic boundary of Q. 

Given appropriate hypothesis on the data, for instance convexity of and twice continuously differ- 
entiable, the classical parabolic theory yields existence and uniqueness for smooth solutions of (|1.5p . 
However classical parabolic theory will only provide estimates involving constants that degenerate as e 
(in the transition from parabolic to degenerate parabolic regime). Our main goal consists in proving stable 
estimates. 

Our first result consists in showing that if the initial/boundary data is sufficiently smooth then the 
solutions of ()1.5|) are Lipschitz up to the boundary uniformly in e > 0. 



Theorem 1.1. (Global gradient bounds) Let G be a Garnot group of step two, Q d G a bounded, open, 
conve^ set and tp € G'^{Q). For 1 > e > denote by Ue G G'^{fl x (0,T)) nC^{Q x (0,T)) the non-negative 
unique solution of the initial value problem (|1.5p . There exists G = G{G, \\f\\c^{fi}) > such that 

(1.6) sup iV^u^l < sup |ViUe| < G. 

nx(0,T) Ox(0,T) 

Remark 1.2. The hypothesis of the Theorem above can be slightly weakened by asking that G is a step two 
Carnot group, tp G G^(f2) but that d^l be required to be convex in the Euclidean sense only in a neighborhood 
of its characteristic locus T,{dn) and that its mean curvature /if^ < ~S < in dfl \ S(9rj). 

Having established Lipschitz bounds, it is easy to show that the right derivatives Xf Uc of the solutions of 
(|1.5p arc themselves solutions of (|3.ip . a divergence form, degenerate parabolic PDE whose weak solutions 
satisfy a Harnack inequality (see [11] and Proposition 13.81 below) . Consequently one obtains G^'" interior 
estimates for the solution which are uniform in e > 0. At this point one rewrites the PDE in (jl.5l) in 
non-divcrgcncc fornix of the total variation flow equation dtu^^ = (VeMe)X'Xjue and invokes the Schaudcr 
estimates in Proposition 14.41 to prove local higher regularity and long time existence. Since all the estimates 
are stable as e — > one immediately obtains smoothness and the consequent global in time existence of the 
solution for the sub-Riemannian case e = 0. 

Theorem 1.3. In the hypothesis of Theorem \l.l\ (or Remark \1.2(l one has that there exists a unique solution 
€ G°°(ri X (0, go)) n i°°((0, oo), G^(f2)) of the initial value problem 



(1.7) 



9tUe = /le m Q = r2 X (0, oo) 
~ p on dpQ 



and that for each k £ N there exists Gk — Gk{G, p, k,Q) > not depending on e such that 
(1-8) \\u,\\cHQ)<Ck. 

Since the estimates are uniform in e and in time, and with respect to e, wc will deduce the following 
corollary: 



■^We say that a set Q G G is convex in the Euchdean sense if exp ^(Sl) C G is convex in the Euclidean space Q. In a group 
of step two this condition is translation invariant. 

^this is possible since for e > the solutions Uf arc sufficiently smooth 
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Corollary 1.4. Under the assumptions of the Theorem \l.l\ as e — > the solutions converge uniformly 
(with all its derivatives) on compact subsets of Q to the unique, smooth solution ua £ C°°{fl x (0, oo)) H 
L°°((0, oo), C^(jl)) of the sub-Riemannian total variation flow (|1.4p in il x (0, oo) with initial data (p. 

Corollary 1.5. Under the assumptions of Theorem ll.ll as T ^ oo the solutions Ug(-,t) converge uniformly 
on compact subsets of D, to the unique solution of the minimal surface equation 

hg: = in D, 

with boundary value ip, while uq = liine_).o Me € C°°{fl)nLip{D,) is the unique solution of the sub-Riemannian 
minimal surfaces equation ho = in fl, with boundary data ip. 



2. Structure stability in the Riemannian limit 

The Carnot-Carathcodory metric is equivalent to a more explicitly defined pseudo-distance function, that 
we will call (improperly) gauge distance, defined as 

(2.1) |xp'^'=^5]|x,|^, anddoix,y) = \y-'x\. 

k=l 1=1 

The ball-box theorem in [33] states that there exists A = A{G, ctq) such that for each x G G, 

A-^\x\ < do{x,0) < A\x\. 

If X e G and r > 0, we will denote by 

B{x,r) = {yeG\ do{x,y) < r} 

the balls in the gauge distance. For each e > wc also define the distance function d^ corresponding to the 
Riemannian metric a^, 

d,{x, y) = inf{ / lyU^ {s)ds with 7 : [0, 1] ^ G 
Jo 

a Lipschitz curve s. t. 7(0) = a;,7(l) = y} 
as well as the pseudo-distance dG,e{x,y) = N^{y~^x) with 

(2.2) N-{x)= -?+nhn|^( ^ 4)^."^ ^ x^X. 

d{i) = l ^ i=2 d(k)=i d(i)>2 ^ 

Note that in the definition of d^^ , if the curve for which the infimum is achieved happens to be horizontal then 
de{x, y) = do{x, y). In general we have sup^^g de{x, y) = do{x, y) and it is well knowrQ that (G, d^) converges 
in the Gromov-Hausdorff sense as e -H- to the sub-Riemannian space (G, do). The ball-box theorem in [39] 
and elementary considerations yield that there exists A = A(G, cto) > independent of e such that for all 
x,y &G 

(2.3) A-'^doAx, y) < de{x, y) < AdG,e{x, y) 
(see for example [TT]V 



'See for instance |30| and references therein 
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2.1. Stability of the homogenous structure as e — > 0. If G is a Carnot group, de is the distance function 
associated to Ce, we will denote 

B,{x,r) = {y€ G\deix,y) < r}. 
If we denote by dx the Lebesgue measure and by |rj| the corresponding measure of a subset f2, then Rca and 
two of the authors have recently proved in [11] that 

Proposition 2.1. There is a constant C independent of e such that for every a; € G and r > Q, 

\B,{xM\<C\B,{x,r)\. 

Having this property the spaces (G, d^, dx) are called homogenous with constant G > independent of e 
(see HO])- 

Let r > and consider the space G = G x (0, t) with its product Lebesgue measure dxdt. In G define 
the pseudo-distance function 

(2.4) d^{{x,t),{y,s)) =max(d,(x,?/),^|t-s|). 

Proposition 12.11 tells us that (G, de, dxdt) is a homogeneous space with constant independent of e > 0. 

In the paper [11] it is also shown that a Poincare inequality holds with a choice of a constant which 
is stable as e — > 0. The stability of the homogenous space structure and of the constant in the Poincare 
inequality arc two of the key factors in the proof of Proposition [3?8| 

2.2. StabiUty in the estimates of the Heat KerneL The results in this section are of independent 
interest and concern uniform Gaussian estimates for the heat kernel associated to elliptic regularizations of 
the Carnot group sub-Laplacians. They will be used in this paper in connection to the Schaudcr estimates 
and the higher regularity of the total variation flow. We will deal with non-divergence form operators similar 
to those in (|1.3p . but with constant coefficients. More precisely we will consider the operator 

n 

(2.5) = dt~Yl a^^'^l 

where A = {aij)ij^i_,,,n is a symmetric, real-valued n x n matrix, such that for some choice of constants 
A, Gi, G2 > and for all ^ e R" one has 

(2.6) + Gi 5: < E a.M, < A E + ^2 E ^f- 

d{i} = l d{i)>l i-j=l d{i) = l d{i)>l 

The point here is that the constants Gi , G2 , A are independent of e and provide A— uniform coercivity of 
(|2.5|) in the horizontal directions. Formally, in the sub-Riemannian limit e — > the equation becomes 

m 

(2.7) LA=dt- a,,Xaj. 

In order to ensure that the operator L^^a (reps. La) is uniformly elliptic (reps, subelliptic) , we will assume 
that the matrix A = (oij ) belongs to a set of the form 

A/a = {A : A is a symmetric n x n, real valued constant matrix, 

satisfying (|2.6p for some choice of Gi and G2} 

for some fixed A > 0. 

Heat kernel estimates in Nilpotent groups are well known (see for instance |33j and references therein). 
We also refer to [5] where a self contained proof is provided. In our work we will need estimates which are 
uniform in the variable e as e — >■ 0, in the same spirit as the results in [TB]. We will denote T^^Aix^t) the 
fundamental solution of (|2.5p . with matrix (oy ) in A/a, and TA{x,t) the fundamental solution of (|2.7p . 
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Proposition 2.2. There exists constants Ca > depending on G,cto,A but independent of e such that for 
each e > 0, X €z G and t > one has 

(2.8) C7^ ^ <T,A(x,t) < Ca ^■ 

For s £ N and k— tuple {ii, . . . ,ik) G {1, . . . , n}^ there exists Cs,k > depending only on fc, s, G,ao,A such 
that 



(2.9) mxi---xir,A)ix,t)\<c,M~'-'^'-^^^^Q^^^^^ 

for all X £ G and t > 0. For any Ai,A2 G Ma, s £ N and k— tuple {ii, ... ,ik) € {1, ■ . . j™}*^ there exists 
Cs,k > depending only on k, s,G,aQ, K such that 

(2.10) \{dtX,,---X,,V,^A,){x,t)-{dtX,,---X,,T,,A,){x,t)\ < 



< \ \A,~A2\\GsM~'~''^^-, 



\BMVt)\' 

where \\A\\^ :-E",=i4- 

Moreover, for s G N and k— tuple {ii, . . . ,ik) G {1, . . . , to}'"" as e — > one has 

(2.11) Xi^ ■ ■ ■ Xi^d^T^^A ■ ■ ■ Xi^dlTA 

uniformly on compact sets and in a dominated way on all G. 



The proof of our result is directly inspired from the proof in |16| . where the time independent case was 
studied by two of us, and where a general procedure was introduced for handling the dependence on e and 
obtain independent estimates in the sub-Laplacian case. 

Let us consider the group G = G x G defined in terms of n new coordinates y G G, so that points of G 
will be denoted x = {x,y). Denote by Yi, . . . , F„ a copy of the vectors Xi,. . . , X„, defined in terms of new 
variables y. The vector fields {Xi) and {Yi) are defined in the product algebra G ~ G x Q. 

In G we will consider two different families of sub-Riemannian structures: 

• A sub-Riemannian structure determined by the choice of horizontal vector fields given by 

(2.12) {Xi, ■ ■ ■ , A"m+n) = (Xi, • . . , Xm, Yi, . . . , Yn). 

The sub-Laplacian/heat operator associated to this structure is 

m-l-n m n 

(2.13) Lo.A = dt~Yl ^* - II a^jX^Xj ~ ^ a^jYiYj. 

i,j = l i,J = l i,J = l 

We complete the horizontal frame to a basis of the whole Lie algebra by adding the commutators 
of the vectors {Xi): {X„i+i, . . .Xn). We denote the exponential coordinates of a point x around a 
point xq in terms of the full frame through the coefficients {v^ , w^) which are defined by 

ni n n 

x = expC^tfiX., + Y,w^Y^+ E «°^0(^o). 

i—1 i—1 i—m-\-l 

• For every e G [0,1) consider a sub-Riemannian structure determined by the choice of horizontal 
vector fields given by 

(2.14) {XI, • • • , Ar^_,_„) = {Xi, . . . , Xm, Yl, . . . , Fm, -|- Yfn+l, . ■ . , X^ + Yn). 
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The sub-Laplacian/heat operator associated to this structure is 

— 1 i J — 1 i>m or j>m 

Analogously, the exponential coordinates associated to L^,a will depend on the given horizontal 
frame (|2.14p and the family (Xm+i, ■ ■ ■ Xn)- The coordinates of a point x are the coefficients w| and 
wf satisfying 

m m n n 

(2.15) x = exp(^^;|X, + ^<y,+ ^ <(X,^+y,)+ E ^^O^*")' 

i—l i—1 i—m+1 i—m+1 

We denote by Vq^a and Fe^a the heat kernels of the corresponding heat operators. In both structures we 
define associated (pseudo)distance functions d[){x, y) and d^{x, y) that are equivalent to those defined in (|2.ip 
and that assign unit weight to the corresponding horizontal vectors while {X„i+i, . . -Xn) will be weighted 
according to their degree d( j). 

Definition 2.3. For every e > and x^xq £ G define 



J,(i,Xo)=El"^^l+El<l+ E min(K|,K|iM^))+ 



i/d(i) 

i—l i—l i— m+1 i—m-\-l 

For e = and x^xq G G define 

n 

jo(i,xo)=5:(k°r/'^«+k°|W) 

We will denote by and Bq the corresponding metric balls. 

Lemma 2.4. There exists constants C\ > depending on G, ctq, A feitt independent of e such that for each 
e > 0, X G G and t > one has 



g-L-A 1 _ g Ca* 

(2.16) ^ < Fe A(S,i) < Ca^ 

i^or s e N and k— tuple {ii, . . . ,ik) G {1, . . . , n + m}'^' i/iere exists Cs,k > depending only on k, s, G, ag, A 
such that 

(2.17) |(a..x-...x.r„.)(..,)l<ft,r.-'/^^;^ 

/or all X E G and t > 0. Moreover, for s £ N and k— tuple (ii, . . . , ik) € {1, . . . , to}'"' as e — > one /las 

(2.18) A^, • • • Xldt^^,A -> AO . . . a,^f o,A 
uniformly on compact sets, in a dominated way on all G. 

Proof. In order to estimate the fundamental solution of the operators L^^a in terms of Lq^a, we define 
a volume preserving change of variables on the Lie algebra and on G: Define : G ^ G a.s F^{x) = 
exp(Te(log(i)) with 

T,{X°) = X^ ioi i = l,...,m + n. 
By definition the relation between the distances do and is expressed by the formula 

(2.19) d,{x,xo) = do{F,{x),F,{xo)). 
Analogously we also have 

(2.20) f,,A(x,t) =fo,A(i^.(x),t), 
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X,r,,A{x,t) = X^fo.A{Fe{x),t), for z = I,-- - ,m 

and 

{Xt +Y^T,^A{x,t) =Y,foMFe{x),t), iori = m + l,--- ,n, 

with similar identities holding for the iterated derivatives. 

In view of the latter, assertions (|2.16p and (|2.17p immediately follow from the well known estimates of f 
(see for instance there references cited above or [5]). The pointwise convergence (|2.18p is also an immediate 
consequence of (|2.19|) and (|2.20p . In order to prove the dominated convergence result we need to relate the 
distances do and d^'. Expressing the exponential coordinates w|, in terms of vf, one easily obtains 

m n 
2—1 2— m+1 

SO that for al0 x,xo £ G 

(2.21) do{x,XQ) - Co < d^{x,xo) < do{x,xo) + Gq 

where Co is independent of e. The latter and (|2.17l) imply dominated convergence on G. □ 

Lemma 2.5. For s G N and k— tuple {ii, . . . ,ik) G {l,...,m}'^ there exists Gg.k > depending only on 
k,s,G,ao,A such that depending on G,ao,A but independent of e such that for each Ai,A2 G A/a, e > 0, 
X € G and t > one has 

(2.22) lid^X^^ ■ ■■Xlf,,A,)ix.t) - [dtXl ■ ■■Xlf,,A.){x,t)\ < 

Proof. In view of (|2.19p and (|2.20p it is sufficient to establish the result for f q.Ai and f 0,^2 1 thus eliminating 
the dependence on e. 

Although the vector fields X^, i ~ 1, ...,m + n are not free, we can invoke the Rothschild and Stein 
lifting theorem [42] and lift them to a new family of free vector fields in a free Carnot group. For the sake 
of simplicity we will continue to use the same notation Xf, i = 1, ...,m + n to denote such family of free 
vector fields. To recover the desired estimate (|2.22p for the original vector fields one needs to argue through 
projections down to the original space, exactly as done in [39| and [16] . A standard argument (see for instance 
[7]) yields that for every A G Ma, if we set A as in (|2.13p then there exists a Lie group authomorphism Fa 
such that 

( {A'/%,X^){u o Fa) = {X,u) o Fa, 

and 

foA^,t) = |detAi/2|fo,7(FA(S),<) 
where / denotes the identity matrix. The authomorphism Fa is defined by 

m+n 

Fa{x) = exp(rA(log(i)), with Ta{X^) = ^ for i = 1, . . . ,m 

and is extended to the whole Lie algebra as morphism (see the Appendix for more details). As in (|2.15p . we 
denote by {vi, Wi) the canonical coordinates of Fa^ {x) around Fa2 (x), then the Mean Value Theorem yields 

{d^xl ■ ■ ■ xlrj){FA, {x), t) - {dtxl ■ ■ ■ xlrj){FA, (x), t) = 



''This estimate indicates the well known fact that at large scale the Ricmannian approximating distances are equivalent to 
the sub-Riemannian distance 



10 L. CAPOGNA, G. CITTI, AND M. MANFREDINI 

n 

i=l 

By the result in the Appendix, the operators ViXi+WiYi arc zero order differential operators whose coefiicicnts 
can be estimated by ||Ai — yl2||. The conclusion follows by virtue of Proposition 15 .11 □ 



We conclude this section with the proof of the main result Proposition [ 

Proof of Proposition \2.3[ From the definition of fundamental solution we have that 

ro,A(a;,i)= / ro,A{{x,y),t)dy, and r^^A{x,t) = / T^^A{{x,y),t)dy, 
Jg Jg 

for any x G G and t > 0. In view of the (global) dominated convergence of the derivatives of T^A to the 
corresponding derivatives of f o,a as e — ?> 0, we deduce that 

^e.A{{x,y),t)dy ^ / To^A{{x,y),t)dy 
G Jg 

as e — )■ 0. The Gaussian estimates of ^e,A follow from the corresponding estimates on f and the fact that 
in view of 

de{{x,y), {xQ,yo)) > da{{x,y), (a;o,2/o)) - Co > do(x,xo) + do{y,yo) - Co > 

> de{x, Xo) + de(y, yo) ~ 3Cq 

Indeed the latter shows that there exists a constant C > depending only on G, (Jq such that for every 

X eG, 

e * dy < Ce * e * dy < Ce * 

JG JG 

The conclusion follows at once 



3. Gradient estimates 

In this section we prove Theorcm ll.il The proof is carried out in two steps: First we use the maximum 
principle to establish interior bounds for the full gradient of the solution Viw of (|1.5p with respect to 
the Lipschitz norm of u on the parabolic boundary. Next, we construct appropriate barriers and invoke the 
comparison principle established in [5] to prove boundary gradient estimates. The combination of the two 
will yield the uniform global Lipschitz bounds. 

3.1. Interior gradient estimates. Recalling that the right invariant vector fields XJ commute with the 
left invariant frame X^, i = 1, . . . , n it is easy to show through a direct computation the following result. 

Lemma 3.1. Let G C^{Q) be a solution to p.3p and denote vq ~ dtu, Vi ~ X^u for i ~ i, . . . ,n. Then 
for every h — 0, . . . ,n one has that Vh is a solution of 



(3.1) dtVh = X^{a,jXjVh) = a'r^{V,u,)X^X]vh + d^.a^^ ,u)XIX'^u,XIvh. 

where 

Note that in order to prove L°° bounds on the horizontal gradient of solutions of (|1.5p one cannot invoke 
Lemma 13.11 with differentiation the horizontal frame, because these vector fields do not commute. The 
argument below bypasses this difficulty by using the algebraic relation between right and left invariant 
vector fields. 
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Proposition 3.2. Letu^ e C^[Q) he a solution to ()1.5p with bounded. There exists C = C{G, \\'f\\c^{n)) > 
such that for every compact subset K dd one has 

sup |Viu| < sup(|Viu| + Ic'twl), 

Kx[Q,T) dpQ 

where Vi is the full ai — Riemannian gradient. 

Proof. The proof is similar to the argument in [8l Proposition 5.1] and is based on the trivial observation 
that left-invariant vector fields at the layer k can be written as linear combination of right-invariant vector 
fields at level k and left (or right) invariant vector fields at levels fe -I- 1, fe -I- 2, . . . , r. The first step consists 
in differentiating the equation in (II. 5p in the direction of the highest layer (which commutes with all the 
Lie algebra) and invoke the weak maximum principle applied to p.ip to show 

(3.2) sup V \X,u\ < sup(|Viu| + \dtu\). 

Next one proceeds in a similar fashion, iterating the argument by differentiating the equation in (|1.5p along 
the right invariant vector fields in the layer V^~^. The maximum principle and p.ip yield L°° bounds on 
XI u with d{Xl) = r — 1. Coupling such bounds with (|3.2p one obtains 

(3.3) sup V \X,u\ < sup(|Viu| + \dtu\). 

^x[0-^)d(.o=r-i ^^'^ 

The proof now follows by induction on the step of the derivatives. □ 



3.2. Linear barrier functions. In [8l Section 4.2] it is shown that in a step two Carnot group coordinate 
planes (i.e. images under the exponential of level sets of the form Xk ~ 0) solve the minimal surface equation 
ho = 0. In the same paper it is also shown that this may fail for step three or higher. In the construction of 
the barrier function we will need the following slight refinement of this result, 

Lemma 3.3. Let G be a step two Carnot group. // / : G — > K is linear (in exponential coordinates) then 
for every e > 0, the matrix with entries XfX^f is anti- symmetric, in particular every level set of f satisfies 
K = 0. 

Proof. We need to show that for f{x) 

n 

We recall the expression (jl.2p for the vector fields Xi, d{i) = 1 in terms of exponential coordinates Xi = 
5a;i+X)d(j)=i d(h)=2 ^%^j^xh - Thc Campbell-Hausdorff formula implies the anti-symmetry relation c^^- — — c^^- 
It is immediate to observe that, if d{k) = 2 one has 

(3.5) X,{xk) = ^ dtjXj, and XiXj{xk) ^ c^i, for d{i) = d{j) ^ 1, 

dU) = l 

if either d{i) — 2 ot d(j) = 2 it is easy to check X^Xj(xk) ~ for all fc = 1, . . . , n. Since ^(5^ — r^q^^p-^ 
is symmetric in i,j it follows that 

(3.6) h^^yakyi 5,, - 1 XtXUxu) = 0. 





□ 
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3.3. Boundary gradient estimates. We say that a set Q C G is convex in tlic Euclidean sense if 
exp^^(fi) C is convex in the Euclidean space Q. In a group of step two this condition is translation 
invariant. 

Proposition 3.4. Let G be a Carnot group of step two, fl G G a bounded, open, convex (in the Euclidean 
sense) set and ip € C^(j7). Fore > denote by G C^(Ox (0, T)) H C^(r2 x (0,T)) the non-negative unique 
solution of the initial value problem (|1.5p . There exists C ~ C{G, \ \'f\\c^ (p)) > such that 

(3.7) sup \VeUe\ < sup iViu^l < C. 

dnx(0,T) dnx{0,T) 

We start by recalling an immediate consequence of the proof of [HI Theorem 3.3]. 

Lemma 3.5. For each e > 0, ifUc is a bounded subsolution and is a bounded supersolution of p. 51) then 
Ueix,t) < Veix,t) for all {x,t) e Q. 

Let Me S G^{Q) be a solution of (jl.Sp . and express the evolution PDE in non-divergence form 

n 

(3.8) dtUe =he= alj{\7eu)X^X]u,. 

Set — Ue — so that fe solves the homogenous 'boundary' value problem 



(3.9) 



dtv, = {V,v, + V^ip)XlXp, + b' in = X (0, T) 
~ on dpQ, 



with b'^{x) = alj{\/eVc{x) + \7i(p{x))X^Xjip{x). We define our (weakly) parabolic operator for which the 
function is a solution 

(3.10) Qiv) = atj{V,v, + V,ip)XtX]v, + b'- dtv. 

In the following we construct for each point po = {xq, to) G dfl x (0, T) a barrier function for Q, v^: i.e., 

Lemma 3.6. Let G be a Carnot group of step two and Q G G convex in the Euclidean sense. For each point 
Po = {xo,to) € dfl X (0,T) one can construct a positive function w G C^{Q) such that 

(3.11) Q(w) < in V n Q with V a parabolic neighborhood of po, 

w(po) = and w > in dpV D Q. 

Proof. In the hypothesis that Q is convex in the Euclidean sense we have that every xo S dfl supports a 
tangent hyperplane P defined by an equation of the form n(a;) = X]r=i ~ ^ with 11 > in fi, n(a;o) = 0, 
and normalized as X]!i(i)=i ~ ^- Following the standard argument (sec for instance [34l Chapter 10]) we 
select the barrier at (a;o,to) G 9il, x (0,T) independent of time with 

(3.12) w = $(n) 
with $ solution of 

(3.13) $" + i/($')^ = 0, 
in particular 

(3.14) <^>{s) = -log{l + ks), 

V 

with k and v chosen appropriately so that conditions p. lip will hold. We choose a neighborhood V = 
O x (0, T) such that Pn O n dVl = {xq}. By an appropriate choice of k sufhciently large we can easily obtain 
Wfi{po) = and > Ve in dpV n Q. 
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To verify Q{'We) < we begin by observing that w satisfies 

(3.15) Q{w) = ^'at,{VeW + V,ip)X^X;n+^T + K, 

witli J" = ajjiVeW + V^(p)XfwX^w. 
We will show: 

(Claim 1) atjiVeW + \',ip)X^X^n < 

(Claim 2) _|^^ + 5, <0 

in a parabolic neighborhood of pq. 

The first claim holds with an equality as af^- is symmetric and X^XjH is anti-symmetric in view of Lemma 
13.31 To establish (jClaim 2p we first note that Lemma [3.31 implies 

^<max( ^ ale' ^ 4) < \V = 

d(i) = l d(fe) = l 

= E (°' + E E al<C{G){l + e% 

d{i) = l ^ d(fc)=2,d(j) = l ^ d{k)=2 

for some constant C{G) > 0. Consequently, for $' >> 1 sufficiently large one finds 

(3-16) :F>- .^'^^""j! >C{G)- — — - > C(G)e2 > 0, 

with C(G) > a constant depending only on G (not always the same along the chain of inequalities). In 
view of the definition of be and (|3.13p with an appropriate choice oi f = v{G, e^cj)) > and k = fc(G, 0) >> 1 
in (|3.14p . we conclude 

□ 

In view of Lemma 13.51 a comparison with the barrier constructed above yields that 
(3.18) < < ^(^'^) < c(fc, v), 

inVnQ, with dista-i{x,xo) being the distance between x and xq in the Riemannian metric cri, concluding 
the proof of the boundary gradient estimates. 



The proof of Theorem 11.11 now follows immediately from Proposition 13.21 and Proposition 13.41 

Having established uniform global Lipschitz bounds one now notes that equation p.ip satisfies horizontal 
coercivity conditions uniformly in e > 0. Such conditions are among the main hypothesis of the Harnack 
inequality in |llj . In this paper, G. Rca and two of the authors have proved that given a homogenous 
structure and a Poincarc inequality, then a sub-elliptic analogue of Aronson and Serrin's Harnack inequality 
[3] for quasilinear parabolic equations holds. As a consequence one obtains for some a € (0,1) that the 
solutions Ue to (jl.Sp satisfy G^'" Holder estimates, uniform in e € (0, 1), 
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Definition 3.7. Let < a < 1, Q C and u be defined on Q. We say tliat u £ C^xiQ) if tliere exists 

a positive constant M sucli tliat for every {x,t), (a;o,io) G Q 

(3.19) \uix,t) -u{xo,to)\ <Md'^{{x,t),{xoM))- 

We put 

|u(x,i) - w(a;o,<o)| , , , 
||u||c° (Q) = sup ^= hsup|w|. 

(x,t)^[xoM) df{{x,t),{xo,to)) Q 



Iterating this definition, if A; > 1 we say that u E C'^'^{Q) if for all z = 1, . . . , m <E C'^ j^'" {Q). Where we 
have set C°;^(Q) = Cf^xiQ)- 

Corollary 3.8. (Interior C]^" estimates) In the hypothesis of the previous results, letting K be a compact 
set K CC Q, there exist constants a G (0, 1) and C = C{K, a) > such that for all i = 1, . . . , n one has 
that V = Xfu satisfies 

M\cf^^{K) + \\'^eV\\mK)<C, 

uniformly in e E (0, 1). 

4. Regularity properties in the C*-''" spaces 

In this section we will prove uniform estimates for solution of (|1.3p in the C'^'x Holder spaces. This 
is accomplished by using the uniform Gaussian bounds established in Section 12.21 to develop new uniform 
Schauder estimates for solutions of second order sub-elliptic differential equations in non divergence form 

n 

L,u = dtu~ al^{x, t)X'^X]u = 0, 

i,j=l 

in a cylinder Q = H. x {0,T) that are stable as e — > 0. As usual we will make use of the associated linear, 
constant coefhcient frozen operator: 

n 

where {xq, to) S Q. 

As a direct consequence of the definition of fundamental solution one has the following representation 
formula 

Lemma 4.1. Let w be a smooth solution to L^w — f in Q C M"+^. For every cj) G C^{Q), 
(4.1) {wcl)){x,t) - / r^^^ ,^)((x,O,(y,r))(L,(,„,t„)-L,)(u;0)(2/,r)dydT+ 

n 

^x„,to)ii^, t), {y, r)){fc^ + wL,cf> + 2 ^ {y, r)X>X|0) {y, T)dydT, 

where we have denoted by F^^^ the heat kernel for of L^ j^xa.to)- 

We explicitly note that for e > fixed the operator L^ {^xo.to) is uniformly parabolic, ils heat kernel can 
be studied through standard singular integrals theory in the corresponding Riemannian balls. Hence, as 
noted in one is allowed to differentiate twice the kernels defined in (|4.1[) with respect to any right or 
left invariant vector field. 
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Proposition 4.2. Let < a < 1 and w be a smooth solution of L^w = / G C"xiQ) ''^ cylinder Q. Let 
K he a compact sets such that K <Z<Z Q, set 26 ~ do(K, dpQ) and denote by Ks the 5— tubular neighborhood 
of K . Assume that there exists a constant C > such that for every e e (0, 1) 

There exists a constant Ci > depending on S, a, C and the constants in Proposition \2.S\ such that 

\\w\\c2.o,^^^ < C\ (||/||c°^-(if,) + \Mcl;^iKs)) ■ 

Proof. The proof follows the outline of the standard case, as in Theorem 4, Chapter 3], and rests crucially 
on the Gaussian estimates proved in Proposition 12.21 and the fact that the functions [L^j^^^ to) ~ ^e) {w(t>), 
and (/(/) + wL^cj) + 2aljXfwXj(j)) are Holder continuous. 

Choose a parabolic spher^ Bc,s CC K where i5 > will be fixed later and a cut-off function ip G C^(M"+^) 
identically 1 on s/2 and compactly supported in Be,s. This clearly implies that for some constant C > 
depending only on G and ao, |Ve^| < CS~^, \L'^(j>\ < C5~^, in Q. Next, invoking (|4.ip one has that for 
every multi-index / = (ii, 12) S {1, . . . ,m}^ and for every (a;o,io) G B^^g 



(4.2) Xt^X^w^){xo,to) = / Xt^XtTl^^,^^i; (y, r))|(,„,,„) - L,) {w <j))iy,T)dydT+ 

JQ 



I^XlXt^Tl^^^,J;{y,r)) ^/0 + u,L,0 + 2 |] at^X^wX^^c^ {y,T)dydT. 



The uniform Holder continuity of a^^ , with Proposition 12.21 and Lemma 12.51 yield 



<Cd'^{{x,t),ixo,to)) 



|Se(0,V^^)| 

with C > independent of e. In view of the latter, using basic singular integral properties (see [26] ) and 
proceeding as in j291 Theorem 2, Chapter 4], we obtain 



(4.3) 



KX^J'lo,to) (■' - L,M)){w 0)(y, T)dydT 



<Cl (L, -Le,(:ro,to))(^^</') ^„ ,^ , 
i-.j 

< Ci5"\\a: 



„llc°^(S..,)l|W0llc^^-(B,,,). 

where Ci, and Ci are stable as e — >■ 0. Similarly, if (f> is fixed, the Holder norm of the second term in the 
representation formula ()4.2p is bounded by 



(4.4) 



KK^xo.to)((^o, to), {y, r)) (/0(y, r) + i«L0(y, r) -f 2at^XtwX^cl>)dydi 



< C2 



c 



''That is a sphere in the group G = G X R in the pseudo- metric defined in 112.41 1. 
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From (i42|) . (|43l) and (j^i]) we deduce that 



Choosing 6 sufficiently small we prove the assertion on the fixed sphere B^ s The conclusion follows from a 
standard covering argument. 

□ 

A direct calculation shows the following 

Lemma 4.3. Let k ^ N and consider a k~tuple (ii, . . . ,ik) G {!,..., m}*^. There exists a differential 
operator B of order k + 1, depending on horizontal derivatives of ajj of order at most k, such that 

n 

XI ■ • - L,) = J2 {<j - at,{xo,t,))xi ■ ■ ■ X^X^X^ + B. 

Proposition 4.4. Let w be a smooth solution of L'^w ^ f on Q. Let K be a compact sets such that K dd Q, 
set 25 = dQ{K,dpQ) and denote by Kg the S^tubular neighborhood of K .Assume that there exists a constant 
C > such that 

\\<^i3\\c^-^{Ks) - 

for any e S (0, 1). There exists a constant Ci > depending on a, C , 5, and the constants in Proposition 
\2.Si but independent of e, such that 

Proof. The proof is similar to the previous one for k ~ 1. We note that differentiating the representation 
formula (|4.ip along any fc— tuple {ii, . . . ,ik) G {1, ■ • ■ , m}'', and using Lemma l473l yields 

Xt^ ■ ■■Xl{wc^){x,t) = I Tl^^^,^){{x,t), {y,T))(al^ - <^.(xo,to))n • ■ ■ X^XtX^^iw <ly){y,T)dydT 



Q 



^.oA,)((^'t)'iy'^))Bi^^)iy,r)dyd 



T 



/ yi,,to){{^,t),{y,T))Xl - ■ ■ Xl(fcj,{y,T) + wL,^{y,T) + at^X'^wX'^^^dydr, 



where (j> is as in the proof of Proposition 14.21 and i? is a differential operator of order A: + 1 . The conclusion 
follows by further differentiating the previous representation formula along two horizontal directions Xj_^ Xj^ 
and arguing as in the proof of Proposition 14.21 □ 



Proof of Theorem ] We will prove by induction that for every k Cz N and for every compact set K <Z<Z Q 
there exists a positive constant C such that 

(4.5) < C, 

for every e > 0. The assertion is true if k = 2, by CoroUarv 13.81 and Proposition If (|4.5p is true for k + 1, 
then the coefficients in p.3p satisfy af^ G C^x uniformly as e G (0, 1) and (|4.5p thus holds at the level k + 2 
by virtue of Proposition 14.41 □ 
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5. Appendix 

Let us assume that G be a free nilpotent Lie group, and let Q ~ (B ... (B its associated Lie algebra. 
Denote by Xi, . . . , X„i a basis of the first layer of the Lie algebra, and by X^+i, ■ • • , Xn a list of vectors 
which complete a basis of the tangent space. If u(x) is an homogeneous polynomial of order p and X is a 
differential operator of degree q we will call u{x)X differential operator of degree q ~ p. If aij = {A)ij is a 
real valued, constant coefHcient, m x m positive definite matrix , one can define a Lie algebra automorphism 

m 

Ta:V^^V\ with TA{X^) = 
for every i = 1 , . . . , m and extend it to the whole Lie algebra as a morphism 

m 

TA[Xi, Xj] = [TA{Xi),TA(Xj)] = ^ aihajk[Xh,Xk]. 

hM=l 

Note that Ta can be represented as a block matrix A ~ Ai Q) A2 ■ ■ ■ (S) Ar, with Ai = A and where the block 
Aj act on vectors of degree j, and its coefficients are polynomial of order j of the coefficients of Ai. In 
particular Ta{X) and X have the same degree. Via the exponential map Ta induces a group automorphism 
on the whole group 

Fa:G^G, defined by Fa{x) = exp(rA(log(a;)). 
In terms of exponential coordinates ( around the origin) 

n 

(5.1) a; = exp(^t;,X,)(0) 
one has 

n r 

Fa{x) = exp(^ v,Ta{X,)){{)) = exp(^ ^ «,(^),,Xj)(0), 

i=l s=l d(i)=d(j) = s 

where d{i) denotes the homogenous degree defined as in (|l.ip . Since ViXi is a zero order homogeneous 
operator then also ViTA{Xi) is a differential operator of order zero. 

Proposition 5.1. Let (wi)i=i^...,n denote the canonical coordinates of Fa{x) around Fb{x), and (f j)j=i,...,n 
the coordinates of x around 0, as in (|5.ip . There exists M G N depending only on the group structure; 
constants ci,...,ca/ depending only on the group structure; zero-order differential operators Yi,...,Ya/ 
depending only on the group structure, on the coefficients {vj)j=i „ and their derivatives along vector fields 
up to order r — 1 (but independent on A and B ), and a constant C ~ C(||^||,||i3||,G) > such that 

n M 

(5.2) ^WiX^^^ciYi and\ci\<C\\A- B\\ forl = \,...,M. 
1=1 1=1 

Proof. By definition 

n n n 

^ WjX, = ^ u,Ta{X,) * (- ^ v^Tb{X,)) 

j=l i=l i=l 

where * is the Baker-Campbell-Hausdorff operation (see for instance [H] and references therein). This 
formula allows to represent the left-hand side as a finitcH sum of terms involving commutators (up to order 



'since the group is nilpotent 
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r) of the two terms on the right-hand side. Let us consider separately the terms in this representation 
starting with those involving no commutators, i.e. 



1=1 1=1 s=l d(i)=d{j)=s 



Clearly this expression is in the same form as (j5.2[) . Next wc consider the second term in the Bakcr-Campbell- 
HausdorfF sum, i.e. that involving commutators of the form 



i=l 



h=l 



p,s=l d{i)=d(j)=d(h)=d{k)=s 



VtXj,VhXk 



-(^{As)ij{Bp)hk - iAp)hkiBs)ij^ ViXj,VhXk 



E E 

P,s=l d{i)=d(j)=d(h)=d{k)=s 

It is evident that the coefficients of the commutators above are Lipschitz functions of A and B, vanishing 
for A = B, and consequently the whole sum satisfies (|5.2p . To conclude the proof we note that all further 
terms of the Baker-Campbell-Hausdorff sum involve operators as in (|5.2p and zero order operators of the 
form X]"=i '^iTA{Xi) or — X]"=i '^i^'B^Xi). Such commutators give rise to differential operators of zero order 
whose coefficients are polynomials in A and B and vanish when A = B. □ 
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